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The problem of inferring the outcome of a simultaneous measurement of two non-commuting 
observables is addressed. We show that for certain pairs with dense spectra, precise inferences of 
the measurement outcomes are possible in pre- and postselected ensembles, and if the selections 
involve entangled states with some other system. We show that the problem is related to the 
problem of assigning weak values to a continuous family of operators, and give explicit examples 
where this problem is solvable. Some foundational implications are briefly discussed. 



A joint measurement of two non-commuting observ- 
ables can be understood as the simultaneous coupling of 
the measured system to two independent instruments, 
each of which is a good probe of either observable when 
coupled individually Independently of whatever 

meaning-if any-one may choose to attach to the instru- 
ment readings in this context, no state preparation can 
be realized so as to systematically ensure a definite out- 
come in any subsequent measurement of this type 
Instead, the mutual back-action of the instruments en- 
sures an uncertainty relation for the joint outcomes with 
a lower bound that is twice that given by the standard 
uncertainty principle Q . 

Fully accepting these facts, consider a scenario in the 
same spirit of the so called "King's problem" 0, H, @, 0, 11] : 
Alice sends a particle to Bob, which he then subjects to a 
joint measurement of a non-commuting observable pair- 
say X and p ; finally, he returns the particle to Alice. 
Suppose Bob chooses not to reveal his readings. The 
question is: are there conditions under which the readings 
can still be inferred, precisely, by Alice? 

In this letter we show that this question can be an- 
swered in the affirmative for certain operator pairs, and 
if she makes appropriate measurements both before and 
after the particle has passed through Bob's instruments. 
Moreover, she will need entanglement resources in order 
to achieve the task, not unlike other quantum inference 
problems where entanglement is required for optimality 
(e.g., [1, Therefore, we beheve the results point to 
an interesting connection between quantum uncertainty, 
non-locality and ultimately, classicality. These questions 
will be addressed briefly at the conclusion. For the mo- 
ment, we begin with the formal statement of the problem: 

Joint Measurement Inference Problem: Let (Ai, A2) 
be a pair of observables of a system with no common 
eigenstates, and U{qi, ^2) a unitary operator of the form 

^^(91,92) = exp[i(iigi +1292)] (1) 

{h = 1). When (91,92) are operator- valued, U{qi,q2) is 
the unitary evolution operator describing an impulsive 



simultaneous interaction of the system with two separate 
instruments (/i,/2) described by canonical pairs (gi,7ri) 
and (92,7r2)- Given this interaction, the problem is to 
find a realizable set of conditions not involving the in- 
struments such that for an arbitrary initial state p of 
(/i,/2), the conditional post-interaction density matrix 
p' of the instruments, is related to p through a quantum 
operation with a unitary normalized Kraus operator 



FpF\ F = exp [i(Q!igi + (1292)] 



(2) 



for some p-independent pair of real numbers (011,02). If 
the problem has a solution, then the conditional probabil- 
ity distribution of the "pointer variables" (tti, 7r2) satisfies 
'^(I'l, T^2\p') — ^{'Ki — oi, 7r2 — a2|/5) , and the conditional 
outcomes can be ascertained to within the uncertainties 
in (7ri,7r2). Thus, the outcomes (Q!i,a2) can be unequiv- 
ocally ascertained for any single trial in the "sharp" limit 
(Atti, A7r2) 0. 

We discuss solutions to the inference problem involv- 
ing complete pre- and postselection measurements, in 
which case the Kraus operator is F oc ('0/|J7((7i, 92 )!''/'») 
with pure initial and final states ji/'i) and IV'/)- Re- 
calling the definition of the weak value of a quantum 
mechanical observable 11 , l2|, Ayj = {tpf\A\ipi) / {'ipf\tpi), 



the inference problem is then equivalent to finding a pair 
(iV'i), IV'/)) such that weak values e'*^"!'^''""^'^^) can be as- 
signed to all elements of the continuous set of operators 
{C/((7i, 92)1(91, 92) S M?}. Alice's task should then be 
to perform pre- and postselection measurements, yield- 
ing initial and final states for which this assignment is 
possible. 

Let us then consider the case of the canonical variable 
pair {x,p) for a particle in one dimension. A natural first 
guess at a solution would seem to be a pre-and posts- 
elected ensemble defined by initial and final eigenstates 
of p and X respectively. However, applying the Baker- 
Hausdorff lemma to the exponential in (p|e*^^'^^~'"'^'^^''|a;), 
we find the Kraus operator F = g*2;qi+ip(j2 ^^qiq-i ^ differ- 
ing from the desired form by the additional term e'^^''^. 
This term represents the back-action between the two 



2 



measurements, as it generates the canonical transforma- 
tion (Jtti, (57r2) = 5(92,91); in the absence of correlations 
in the preparation of the instruments, this back-action 
term can be shown to enforce the uncertainty relations 
in the final outcomes (ttJ^, ttj) 



AvriATTj > ATTiATTa 



1 



16A7riA7r2 



> 1/2. 



(3) 



The conditional uncertainty bounds now coincide with 
standard uncertainty relations, suggesting that even with 
additional information from postselection, the standard 
limit is still unbreachable due to inevitable back-action. 

Note, however, that we have only looked at complete 
measurements on the system. So suppose instead that 
the selections involve measurements of general observ- 
ables of the system and an ancillary canonical system, 
with canonical variables {xa,Pa)- In particular, consider 
the conjugate pairs {x^,p_^) defined by 



i+ = (i + Xa)/2 



P+ =P + Pa, (4a) 

P_^ip-pa)/2, (4b) 



so that X = x^ + x_/2 and p = p^/2 + p_. The linear 
combination x qi+pq2 can then be expressed as the sum 
of the two commuting combinations: {x_ qi +p^q2)/2 and 
+ P-Q2, each involving commuting operators; thus, 
we can perform the factorization 



J{xqi+pq2} 



\ 91 -1-33+ 92 ) J.('2^ 91 92 ) 



(5) 



Taking eigenstates of the two terms in the factorization, 
for instance an initial eigenstate of {x_ , p^ ) and a final 
eigenstate of {x^,p_)^ we obtain the unitary normalized 
Kraus operator 



F 



{x.pje 



i{xqi+pq2 



Jixqi+pq2) 



{x^p_\x_p_^) 



(6) 



with X = x_^_+x_/2 andp = p_ +p_^/2. With the absence 
of the back-action term, the ensemble with \ipf) = \x_p_^) 
and — \x^p_) solves the inference problem for {x,p). 

This solution ensemble reveals another surprising re- 
sult when we consider, with otherwise the same condi- 
tions, the joint measurement of all four canonical vari- 
ables of the system plus the ancilla with four different 
instruments. The corresponding Kraus operator then 
involves the quantity |e*(*«i+P«2+*"«3+P"«'') 
Disentangling the exponential, we obtain 



_ „i(2;9i+P92+2:a9.'i+Pa94) „^(9l94+92 93) 



F = e 



(7) 



where Xa — x^ — x_ /2 and Pa = —{p_— P+/2). It 
is therefore the crossed pairs of commuting observables 
{x,pa) and (p, Xa), that generate a back-action term and 
hence uncertainty relations as those of Eq. ^ for initially 
uncorrelated instruments. These results, as the results 
of our earlier failed attempt, are consistent with what 



could be conjectured to be a new type of uncertainty 
relation for a composite system with two canonical de- 
grees of freedom, allowing for any pre- and postselection 
described by states of the composite system only. In 
the absence of any information about correlations in the 
preparation of the measuring instruments, it appears as 
if the uncertainties in any possible inference of the read- 
ings (tt^j , TTp^ , 7r^2 5 '''P2) from a joint measurement of all 
canonical variables are only constrained by the relation 



A< Att' A< AV >1/4. 



(8) 



It is as if entanglement allowed us to arbitrarily deform 
an uncertainty volume in the "inference phase-space", 
provided the total volume of 1/4 is preserved. 

Prior to further speculation, we turn our attention to 
the assignment of weak values associated with the general 
inference problem. The assignment problem can be cast 
in terms of a complex linear map 



(9) 



where W, henceforth termed the weak value operator, 
is realized by some pair {It/ji) , IV'/)) of non-orthogonal 
states in the total Hilbert of the system and an ancilla. 
The general question is then to determine for a given set 
of pairs (i^*), a^*)), an W such that W : A^''> a^^\ 

Consider the finite-dimensional case first. For the 
space of linear operators acting on a Hilbert space Tis of 
dimension d, we introduce a basis ....,Ej2, orthonor- 
mal with respect to the standard hermitian inner prod- 
uct: i.e., {Ei\Ej) = where (i|i3) = Tv{A^B). An 
arbitrary observable A of the system can then be repre- 
sented as a vector a G through A — ^ ■ OiEi = a E 
with Oi = {Ei\A), and where the dot product is euclidean 
{a b = ctibi)- We will make an exception for the weak 
value operator W, which we choose to expand in terms 
of the hermitian conjugate basis: W — w ■ E'^ , in which 
case Wa is the weak value of Ea, and the weak value of 
A can be written a& a ■ w. We reserve the notation / 
for the vector representing the identity operator, where 
li = Tr(£'i)*. One constraint on w is then I ■ w = 1. 

The question of whether w has other constraints in 
,2 

C besides I ■ w — 1 brings us to the significance of 
entanglement. Suppose that one could write W in the 
form W ~ ^(x/ 1x0 ^'^^ normalized \xi) and \xf) in the 
system Hilbert space. It is then easy to see that W sat- 
isfies the relations = W , WW^ = jPA^, and 

W^W = TT^-T^Hvi- These conditions lead to non-trivial 

\{Xf\Xi)\\ 

constraints involving the components of the weak vector 
w, the simplest of which are: 



1 and w ■ w* 



i(x/ix.)r'>i, (10) 



where 7.^ = TT{EiEj). On the other hand, define two 
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vectors in a Hilbert space Hs x Ha, with dim{Ha) > d 

n 

l<i') = El*)^N)- |^) = ^.£;t|$). (11) 
1=1 

Using these as initial and final states for a pre-and posts- 
election and tracing over the ancilla, we obtain a realiza- 
tion of W with w = z/{z ■ I), where no constraints are 
required for z other than 2;-/ 7^ 0. Thus we see that with 
entanglement, it is possible to realize any complex vector 
w that solves the associated linear- algebraic problem of 
satisfying a'*-* = a'*^ ■ w for a vector set {a^*'}, provided 
the problem is solvable and consistent with I ■ w = 1. 

Turning then to the linear algebra problem associated 
with : U{qi,q2) e^("i9i+«292)^ 

we now show that the 
problem has no solution if I3e = Ai cos + A2 sin 6 has a 
discrete spectrum with in any finite subinterval of [0, tt] , 
as is always the case in the finite-dimensional case. First 
note by uniqueness of the exponential expansion, that we 
must have W : Bg ^ (ig where Pg = [ax cos 9 + a2 sm0) 
for all integer powers fc; on the other hand, by the 
Cayley-Hamilton theorem, we also know that Bg annihi- 
lates its characteristic polynomial pg{z) = det(zl — Bg), 
and hence, by linearity of weak values, we conclude that 
PeiPe) = and hence that (3g must be one of the eigen- 
values bg ^ of Bg. But setting 6* = or 7r/2 we see that 
(ai, a2) must be a pair of eigenvalues (ai, 02) of (Ai, A2). 
The conjunction of conditions can be visualized as the 
intersection of the curve (3g = ai cosO -\- 02 sin with all 
zeroes of pg{z) plotted as a function of 6. If Bg has a 
discrete spectrum in some interval, then either the inter- 
section occurs at discrete values of 0, in which case the 
solution fails, or else there must exist a root behaving 
like bg — aicos^ + 02 sin 6* in that interval. But this 
can only occur if Ai and A2 have a common eigenstate, 
in contradiction with our stated assumptions. Conse- 
quently, solutions to the inference problem can only be 
found in the infinite dimensional case and only for op- 
erators Ai and A2 such that the combination Bg has a 
dense spectrum in some band around aicos^ -I- 02 sin 
for some pair of eigenvalues (ai, 61) of (^1,^2). 

Approximate assignments of the form e*"i'?i+"292+o(ij=) 
for some power s can nevertheless be constructed in the 
finite-dimensional case. A generic form for the leading 
correction is obtained when Bg has a minimal polynomial 
mg{z) of degree s (for all 0), and assuming the linear in- 
dependence of all symmetrized operators Si^m generated 
by the expansion (Ait + A2Y = I]f=o v.(k-iy. ^i'k~i for 
all k < s. Then one can assign W : {Aiqi + ^292)*^ ^ 
{aiqi + 0292)'^ for all k < s with arbitrary (ai, 0:2). The 
leading order term in the exponent is then found to be 
—i'^\q\^me{z)/s\, evaluated at z = (ai cos0 -|- 02 sin0) 
(An example is a pair of orthogonal spin directions, in a 
spin j representation, in which case mg{z) = pg{z) and 
s = 2j + l). 



Proceeding with the continuous variable case for a one- 
dimensional canonical system, we generalize previous re- 
sults for two standard basis sets in the Weyl representa- 
tion (see e.g., 13, 3]): the Heisenberg basis and its 



reciprocal basis, where C, and 77 are two-component 
dense indices valued on the standard symplectic plane 
of dimension two. Throughout, we denote the canoni- 
cal observable pair as f) = {x,p) and use the symplectic 
product notation ^ Ai] = — C^Vi)- The Heisenberg 

basis consists of the translation operators — e"''^^ sat- 
isfying T^fiTQ = fl + C , while the A^ operators form with 
the a Fourier-transform pair 



(12) 



The integration measure is defined as = (27r)~^ J 
etc., with the 27r replaced hy h = 27rft when using units; 
note the orthonormality of both bases with respect to 
this measure, i.e., {ri\r]') — 2TTS{ri — 77'), etc. A summary 
of useful algebraic relations is 

f^f^, = e^^if<;+e , A^A^ = 4e^^^r2(^') , (13) 
fcA„ = e'*^A^+^/2, A^f^ =e^^3A^_^/2, (14) 

where 0i = A C, 02 = 2?? A t]' and <f)z = rj txQ. For the 
expansion of an observable A in either basis 



A 



(C)rc= / S(ry)A,, 



(15) 



the expansion functions a(C), and 0(77) are respectively 
known as the Weyl transform and the Weyl symbol of 
A and form a Fourier transform pair; for the identity 
operator, I{Q = 2n5{C,) and /(ry) = 1. The excep- 
tion again is the weak value operator W , which is ex- 
panded in the Hermitian conjugate of the Heisenberg ba- 
sis W = J^w{C)fl, in which case the Weyl transform 
of w{ri) corresponds to the function w{—(). With these 
conventions, W : A ^ J^a{(^)w{() = J^^a{i])w{ri). 

In the infinite-dimensional case, the weak value assign- 
ment associated with the inference problem involves in- 
tegral equations with solutions obtained by inverting po- 
tentially complicated kernels. For the moment, it will 
suffice to show a generic solution for a reduced set of ob- 
servable pairs, leaving open the question of the general 
set of pairs for which the problem is solvable. We consider 
pairs of the form (p,/(a;)) for functions f{x) satisfying 
conditions determined by the generic solution. Letting 
U{ti,t2) = exp[i{f{x)t2 — pti)], the problem is then to 
find WKit>{C) such that for all (ti,i2), 

e^(**^-''*i) =^i^(t|C)u;,0(C), (16) 



where the kernel K{t\C) = Tr [{7(ti,t2)f/ 



can be shown 
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to take the form 

K{t\C)=d{Ci-ti) I dwe'[»("l^i)*^-^="l , (17) 



test functions, the corresponding Weyl symbol is then 



'•;d„/\d. 



with g{u\Ci) = i J^ids f{u — sCi/2) . If f{x) allows the 
exchange / dC,2 jdu ^ / du / dC2 (see [HI), then it is 
easily verified that w^,^ {Q must be the Fourier transform 
with respect to the C2 variable of a J- function at some real 
root W0(Ci) of the equation g{u\C,i) = 4>. Thus, 



Wk,/,(C) 



=»(«0(Ci)C2-«;Ci) 



(18) 



The function / must therefore be such that g{u\C,i) = (j) 
admits a real branch for all Ci for given 0. These con- 
ditions admit solutions for polynomial f{x), but of odd 
degree only. We have yet to find solutions for bounded 
f{x) (one side or both). 

Proceeding with the realization of solution ensembles, 
it proves instructive to revisit the (x,p) inference prob- 
lem from a constructive viewpoint. From the defini- 
tion of w(C) or equivalently from the above results for 
f{x) = X, we have : — > w{() = e*'''^'» for some 
rjg — {x,p). When the problem is stated in the recip- 
rocal basis, this amounts to finding conditions for which 
w{r]) = 27r5(?7 — r?s), or equivalently, W = A^^. We can 
now work in analogy with the discrete case by introduc- 
ing a maximally entangled state |$o) = / dx\x)s \x)a, and 
two derived basis sets for the combined Hilbert space: 

|$^)=f^|<i>o) and |vI/,,) = A,,|<fo), (19) 

with {'i/rjl^c) ~ e"'^''. Tracing out the ancilla from the 
outer product | for two specific states and nor- 

malizing, we find that 



A 



(20) 



which is the desired result with rjs = rjf + Ci/2. One can 
verify that l^c;), I^',,) are indeed \x_,p_^), \x_^_,p_) with 
C = {x_,P+), V = (a;+,p_)- 

One noteworthy aspect of this solution is the delta- 
function Weyl symbol for W. The ensemble therefore 
has the remarkable property that the weak value of any 
observable is the respective Weyl symbol evaluated at 
ij = rjs- Thus, when observables with classical Weyl sym- 
bols are weakly (hence, unsharply) probed, the condi- 
tional effects on the instruments from this ensemble will 
be indistinguishable from those of a classical system with 
definite phase-space localization. 



121) 



(see also 

Any semblance of standard textbook classicality proves 
illusory, of course, in light of the full repertoire of non- 
local observables that could be probed weakly. Adapt- 
ing equation ^ to the present language by introduc- 
ing the time inversion operation (^'^ = (Cii~C2), we 
find for the composite system the Weyl transform 

w^^'^^CC) = e*(''=^C+''»Oeif^(';')'^ where for the an- 
cilla, 77a = {xa,Pa) = {vf ~ C//2)"^- Notc that time in- 



'v'^^^'S^rj - r]s)5{'q' - rja) . Thus, while 



joint weak measurements yield Weyl symbols for system 
and ancilla observables A and B respectively, weak mea- 
surements of AB will differ from the product of the Weyl 
symbols by terms involving "crossed" derivatives to all 
orders (equivalently, the weak value of AB determines 
correlations in the outcomes of the joint weak measure- 
ment of A and B [l6|). The signature of entanglement is 
therefore a joint Weyl symbol that is essentially a non- 
local object in the combined "phase-space" . 

We conclude with some remarks on interpretational is- 
sues. First, it is hard to ignore the fact that for those 
conditions solving the canonical inference problem, the 
conditionally sharp outcomes {x,p) are in consistent cor- 
respondence with numbers representing physical proper- 
ties of the system in other contexts. To wit: we find 
algebraic correspondence with the initial and final state 
labels through the transformations of Eq. (|4a|) . and at 
the level of weak values, with phase-space functions eval- 
uated at {x,p). These correspondences seem to indicate 
a certain inner consistency of the quantum framework 
which, under special circumstances, allows for a statis- 
tically unambiguous operational assignment of values to 
both canonical variables. Finally, it has been suggested 
earlier [13, [3] that it is not one, but rather a pair of 
state vectors (e.g., our initial and final states) that pro- 
vide the complete description of a quantum system at 
a given instant. The present results provide a good in- 
dication of the extent to which the descriptive power of 
the two-vector framework is enhanced by the property 
of entanglement. It appears that through entanglement, 
the framework becomes flexible enough to incorporate, 
among many other possibilities, the classical description 
of a single system, at least at the kinematic level. We 
believe this realization raises interesting questions, par- 
ticularly regrading our perception of classicality in the 
macroscopic domain. 

The author is grateful to P. Mazur and J. Suzuki for 
helpful discussions. 



version arises naturally from T'i°''|$o) = T):t\^o)- For 
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